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1 . Introduction 


There are two fundamental operations in the load measurement 
problem on an aircraft structure. In the first part, we perform 
a system identification by measuring the responses of the strain 
gages mounted on different locations of the structure from a 
series of known applied load at various specified load points on 
the structure during the calibration stage on the ground. In the 
second part, by using some characterization of the system 
obtained in the first part, we can predict the actual equivalent 
load value and location from the gage measurements during a 
flight. Various known successful approaches and results have 
been reported in the past on the load measurement problem [ 1-3 ] . 

There are two fundamental and intuitively equally justifiable 
linear approaches (arbitrarily denoted as Approach 1 and Approach 
2 in Section 2 ) applicable to the load measurement problem. In 
Approach 1, we model the load value matrix L as dependent 
linearly on the influence coefficient value matrix M measured by 
the gages. In Approach 2, we model M as dependent linearly on L. 
In general these matrices are rectangular, thus it is not immedi- 
ately clear that these two approaches are equivalent. Histori- 
cally, all the work in [1-3] were based on that of Approach 1. 

In Section 2, we shall show that these two approaches are indeed 
equivalent in all cases, and can be proved by the use of the 
modern Singular Value Decomposition (SVD) technique. On the 
other hand, if we only use the more conventional and previously 
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used normal equation technique [1-3] (also called the linear 
regression technique) , then the limitation of this analytical 
technique can only show the validity of Approach 1 when the 
number of gages n is less or equal to the number of loads m. In 
addition, by using the normal equation approach, it is only 
possible to handle these problems with n greater or equal to m 
from the Approach 2 point of view. There are several theoreti- 
cal, practical, and computational consequences to these observa- 
tions. 

At the most basic level of understanding, of course, it is 
theoretically important to know the equivalency of these two 
seemingly different approaches that yield the desired result. At 
the practical algorithmic operational level, the inadmissibility 
of having the number of gages n greater than the number of 
applied loads m in the calibration stage in Approach 1 is not 
fatal. However, as we shall show in Section 3, for the multi- 
stage load estimation technique (which can yield extremely accu- 
rate load predictions) , we will always use more gages than the 
number of loads in the prediction stage. Conventional normal 
equation approach (i.e.. Approach 1) is not possible since a 
crucially needed matrix involved in the processing is singular. 

When the data from the gages are quite linearly independent, 
then there is no significant numerical difference between the use 
of the SVD technique or the normal equation technique. However, 
for highly dependent data, there can be significant advantages 
for the SVD technique. Detailed numerical computations based on 
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practical observed gage measurements and load values are neces- 
sary to verify their differences. The crucial point is that in 
all cases , the SVD approach is always computationally more costly 
as well as numerically more stable. For dimensions encountered 
in the load measurement problems, the additional computational 
cost of the SVD approach is not significant to be of concern. 



2. Two Equivalent Approaches to Load Measurement Evaluations 


Now, we consider two possible basic equivalent approaches for 
load measurements. These approaches are denoted as Approach 1 
and Approach 2 . 

2 . 1 Approach 1 - Linear Dependency of Load Values on Gage Val- 
ues 

Consider 

(1) L = [L 1# L 2 ' L 3 ]* 
the m x 3 load matrix, where 

(2) 1*1 = Lg = [S]_, S 2 , • • • , Sjj] # , 

is the shear vector at m load locations (x^, y^) , 

(3) h 2 = L b = [Bi, B 2 ,..., Bjj] ' , 

is the bending moment vector with its i-th component located at 
y^ given by 

(4) Bi = S^i , i=l,...,m, 

and 

(5) L 3 = I-r = [Ti,..., T m ]\ 

is the torque vector with its i-th component located at x^ given 
by 

(6) T i = SiXi, i = l , ... ,m. 

Let the m x n influence coefficient matrix M denote the 
response of the n gages to the m loads in the calibration pro- 
cess. Specifically, let 
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where each u^, i=l,...,n, represents the normalized response of 
the i-th gages to the m loads. Let the n x 3 dependency coeffi- 
cient matrix fa consists of 

( 8 ) fa = [fai, k 2 > £3]' 
where 

(9) fai - M fai, i=l,2,3, 
or in matrix form 

( 10 ) fa = Ufa. 

For i=l, the n x 1 vector fai yields the dependency of fai = 
fa s , the shear vector, to the linear combinations of the influence 
coefficient vectors (Hi,..., u n ) of M in (7). Similarly, for i = 

2 and 3, fa 2 and fa 3 are related to the bending moment vector fa 2 = 
fa B and the torque moment vector fa 3 = L^. respectively. 

In the calibration process, the matrix Jf as well as fai, fa 2 , 
and fa 3 are available. Define the pseudo-inverse of & as a nxm 
matrix from the "normal equation" point of view as 

(11) M + = (M'M) -1 *!'. 

We note, J! + in (11) is defined if and only if m £ n and all columns 

of J5 are linearly independent. In particular, if n > m, then JJ -1 in 

(11) is not defined. 

Then (9) becomes, 

(12) fai = M + fai i=l ,2,3. 

By using the notation of fa in (8) and fa in (1), (12) can be 
written in matrix form as 


(13) h - M + L. 

In the prediction process, we observe one 1 x n dimensional 
gage measurement vector M (corresponding to the first row vector 
of M in (7)). From (10), the predicted 1x3 load vector is 
given by 

(14) £ = [S, B, J] - & b - 5 IT 1 L = 5 (M'M) -1 M'_L. 

A 

The first component of L yields the predicted shear, 

(15) S - 5 M + La., 

the second and third components of L yield the predicted bending 
moment B = Sy and predicted torque T = Sx, 

(16) B = S y = 5 M + L 2 / 

(17) T = S x = i M + i 3 . 

From (15) and (16) , we can solve for y as 

(18) _ M H + L 2 

y * , 

n U + Li 

and x as 

5 M + L 3 

(19) X = . 

i it Li 

Thus, (15) , (18) and (19) represent the predicted equivalent net 
shear, bending moment location, and torque location of the 
applied load that yielded the measured gage vector M using the 
normal equation approach. 

Now consider the use of the SVD technique via Approach 1. 
Consider a general form of the SVD of the matrix M with rank p 
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£ min{m,n) as given by 

(20) M - U M £m 2'm * 

where U M is a m x m orthogonal matrix, X'jj is a n x n orthogonal 
matrix, and is a m x n matrix where the upper left p x p sec- 
tion is a diagonal matrix with positive singular values (S.V.) 
denoted bya 1 ka 2 ^ ••• £ o p > 0, and the remaining section is 
the zeroth matrix. The pseudo-inverse of the matrix M, is a n x 
m matrix denoted by & ++ , and from the SVD point of view is then 
given by 

(21) *f+ + = Y m £ m + U M ', 

where I M + is a n x m matrix where the upper left p x p section is 
a diagonal matrix with reciprocal of the positive singular values 
of M, and the remaining section is the zeroth matrix. Then by 
using (21) in (10) , we have 

(22) M ++ L - b . 

We note, (22) corresponds to (13) in the calibration stage of the 
previously considered normal equation technique. Then in the 
prediction stage, we have 

(23) L = M b = M M ++ L . 

It is most interesting to note, that the predicted load vector in 
(23) based on the SVD technique has the same form as the pre- 
dicted load vector in (14) based on the normal equation tech- 
nique. Indeed, when m ^ n (i.e., the number of loads is greater or 
equal to the number of gages) , and when the gage measurements are 
quite linearly independent, the pseudo-inverse given by iJ ++ in 
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(21) is equal to the pseudo- inverse given by M -1 " in (11) . Thus, in 
those cases, either the conventional normal equation or the SVD 
methods will yield the same predicted load values. Of course, when 
the measurement values are quite linearly dependent, then the SVD 
approach will be better from the numerical stability point of view. 
As mentioned earlier, when n > m, the normal equation method is 
not applicable for Approach 1 since {T* - in (11) is not defined. 
However, the results of (21) -(23) under the SVD method for 
Approach 1 are valid in all cases including n > m. 


2.2 Approach 2 - Linear Dependency of Gage Values on Load Values 
From a physical cause and effect point of view, it makes 
sense that the responses of the first gage to the m loads are 
given by, 


(24) 

Mi = 

r u n i 

u 12 

• 

• 

— 

r Si 
• 

• 

s lYl 

• 

• 

s l x l ‘ 
• 

• 



• 

L u lm J 


- s m 

• 

s mYm 

• 

s m x m 


= [ Li t L 2 , L 3 ] c 3 — L c 3 . 



In (24), we are describing the gage measurement u^ as a 
linear combination of s^Cn + s^y^c^ + s i x i c 13' which depends 
linearly on the shear, bending moment and torque. In general, 
for all n gages, we have, 

(25) if = [u^,..., u n ] = L [Ci,...c n ] 

= L C, 

where the 3 x n dependency matrix C is denoted by, 

(26) C — [c^,...,c n ]. 


- 8 - 



In the calibration process, as before, M and L are available. In 
the prediction process, as before, we have a measured given 
from (25) as 

(27) M = L C. 

In order to solve for the 1x3 predicted load vector 

(28) L = [ I, £, T ], 

we need to use (25) and (27) . First, consider the use of the 
normal equation technique. Let the psudeo- inverse of L be 
denoted by 

(29) L+ = (L'L)" 1 !/. 

By using (29) in (25) , we have 

(30) C = L + M. 

Substituting (30) in (27), we obtain 

(31) i = 1 1 + H . 

Now, in order to solve for i, we need to multiply both sides of 

(31) by M' from the right and try to take inverse. Unfortu- 
nately, since M is a m x n matrix, if m > n (which is often the 
case), then (MM') is singular and then (MM')" 1 does not exist. 
Thus, in this case, it is not possible to use the normal equation 
technique based on Approach 2. When n £ m, then direct solution 
of (31) yields 

(32) £ = M M' m')' 1 L • 

By comparing (32) of Approach 2 (valid for n £ m) to (14) of 
Approach 1 (valid for m £ n) , we see that only when m = n are 
these two approaches yield identical results. 


- 9 - 



Now, consider solving for £ in (25) by using the pseudo- 
inverse of L based on the SVD representation of L. Specifically, 
consider the SVD of L as given by 

(33) L = U L I L V L ', 

where U L is a m x m orthogonal matrix, is a 3 x 3 orthogonal 
matrix, and E^ is a m x 3 matrix of the form of 

°Zil 0 0 

0 o j'Q 0 

(34) E l = o L3 

m • • 

0 0 0 

where the top 3x3 sub-matrix is a diagonal matrix of singular 
values a L1 £ o L2 Z o L3 > 0 and the remaining (m-3) x3 sub-matrix 
is an all zero matrix. Then (25) becomes 

(35) M = U L Z L V L ' £. 

By using (30) and (33) in (35), we have 

(36) V L Z L + U L 'M = C, 

1 /o 1 0 0 .0 

(37) £ l + = 0 l/o 2 0 . 0 

0 0 l/o 3 . 0 

In particular, we note 

(38) I L + £ L = I 3 . 

By using (36) in (27) , we have 

(39) 5 = L Y L I L + Hl' M • 

Now consider the use of the SVD representation of J1 given by (20) 
in (39). Then we have 

(42) = L V L £ l + JZ l ' H m 3£m/ 
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Direct solution of £ in (42) yields 
(43) 1 = O l+ L . 

By comparing (43) to that of (23) , we see that by using the SVD 
technique, both Approaches 1 and 2 yield the same predicted load 
vector, for all cases of m and n (i.e., m £ n or n > m) . But as 
discussed earlier, by using the normal equation technique. 
Approach 1 is applicable only for m £ n, while Approach 2 is 
applicable only for n £ m. Furthermore, the predicted h in (14) 
under Approach 1 and (32) under Approach 2 are equivalent only 
for m = n. 
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3. Multi-Stage Load Estimation Technique (MUSLET) 


Now, we consider a new proposed MUlti-Stage Load EsTimation 
technique (which we shall denote as MUSLET) for a more accurate 
prediction of the load value. As considered in Sections 1 and 2, 
the load matrix L in equation (1) and the gage measurement matrix 
M in (7) are needed to estimate the dependency matrix in 
Approach 1 and the dependency matrix c in Approach 2 . In either 
case, if we actually know more accurately the true value of the 
load location, then we need to use only those load locations that 
are close (or closest) to the true load location in the calibra- 
tion process to yield a more appropriate b or c which in turn can 
yield a more accurate predicted load value and location. An 
analogous argument to this technique is that in numerical analy- 
sis, only the relevant near by x-y points to the desired x loca- 
tion are used to perform a numerical interpolation. The set of 
far away x-y points may introduce more errors rather than provide 
the desired smoothing effects in the interpolation process. Ini- 
tial numerical experiments based on the HWTS Loads Calibration 
Data of 8/29/76 showed that suppose the 10th calibration load is 
removed from the actual system identification calibration pro- 
cess. By using the remaining 17 load conditions and 15 gage 
measurements, the initial predicted bending moment y location was 
given as 80.3 and the predicted torque x location was given as 
60.2. The true locations were 81.3 and 60 respectively. How- 
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ever, from this initial prediction result, we can use only these 
calibration points with similar or close y values to 80.3. In 
this example, we can use only loads labeled as {7;8;9;11} with 
calibration load values situated at y = 81.3. The newly predicted 
y value now becomes 81.29992. Similarly, the closest x values in 
the calibration stage are given by loads labeled as {3;4;5;9;11; 
15; 17; 18}. The newly predicted x value is now given by 59.84374. 
We note the original predicted values have errors of have errors 
of 1.23% in y and 0.34% in x, while the iterated second stage 
prediction has errors of only 0.0001% in y and 0.19% in x. Thus, 
this MUSLET technique appears to be able to yield significant 
improvement in load predictions. Many variations of the above 
simple proposed iteration schemes are possible and will be 
further investigated. 
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4 . Conclusions 


In this research period from January 1987 to March 1988, we 
have performed various basic research on efficient load measure- 
ment estimation techniques for aircraft structure analysis. In 
Section 1, we presented an over view of the load measurement 
problem. In Section 2, we considered two basic equivalent 
approaches to load measurement evaluations. Under Approach 1, 
the load values are modeled as depending linearly on the measured 
values. Under Approach 2, the measured values depend linearly on 
the load values. By using the modern SVD method, we showed that 
under all conditions of the number of loads m and number of gages 
n. Approach 1 is equivalent to Approach 2. By using the conven- 
tional normal equation (or linear regression) approach. Approach 
1 is only valid for m £ n (which is commonly encountered case) , 
while Approach 2 is valid only for n k m. Furthermore, except 
for the case of m = n, the load prediction formulas under the two 
approaches are not equivalent. 

In many practical flight testing situations, we may not be 
able to use as many gages as those used in the calibration pro- 
cess. Thus, there is much interest in finding the most efficient 
set of n gages to be used for predictions. We have also per- 
formed exhaustive tests on various subsets of the available gages 
and associated load measurements. Preliminary investigations 
show this approach also to be meaningful. Another approach 
based on Distributed Load EsTimation (DISLET) performs a series 
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expansion of the fight measured gage data vector in terms of the 
calibration measured gage data vectors. Coefficients of this 
expansion indicates the amount of load distributions during 
flight relative to the m calibration load points. Further 
analytical and computational investigations are under study. 
Another multi-dimensional cluster analysis relevant in the effi- 
cient reduction of gages, also based on the use of SVD technique, 
has also been under consideration and will be reported on in the 
future . 

Finally, in [4], we have published a basic research paper on 
the effective use of singular values in estimation problems for 
data contaminated by Gaussian noise. The results in [4] yielded 
tighter performance bounds as compared to all previously known 
results in this field. 
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